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Two-point correlation functions of cosmic microwave background polarization provide a physically 
independent probe of the surprising suppression of correlations in the cosmic microwave background 
temperature anisotropies at large angular scales. We investigate correlation functions constructed 
from both the Q and U Stokes parameters and from the E and B polarization components. The 
dominant contribution to these correlation functions comes from local physical effects at the last 
scattering surface or from the epoch of reionization at high redshift, so all should be suppressed if the 
temperature suppression is due to an underlying lack of correlations in the cosmological metric per¬ 
turbations larger than a given scale. We evaluate the correlation functions for the standard ACDM 
cosmology constrained by the observed temperature correlation function, and compute statistics 
characterizing their suppression on large angular scales. Future full-sky polarization maps with 
minimal systematic errors on large angular scales will provide strong tests of whether the observed 
temperature correlation function is a statistical fluke or reflects a fundamental shortcoming of the 
standard cosmological model. 


I. INTRODUCTION 

Two seasons of observational data from the Planck 
satellite have given us the most precise measurement of 
temperature fluctuations in the Cosmic Microwave Back¬ 
ground on the full sky to date These observations 

appear to fit well within the standard picture of our Uni¬ 
verse - Lambda Cold Dark Matter (ACDM). It did, how¬ 
ever confirm several anomalous features in the tempera¬ 
ture fluctuations [4] , which had first been hinted at with 
the COBE-DMR satellite [5] and were later highlighted 
in the WMAP data releases [6]. These anomalies exist 
overwhelmingly at the largest scales of the temperature 
power spectrum, Cj^, with several interesting features 
appearing at multipoles < 30. One feature, the lack of 
two-point correlation at angular separations of 60° and 
above, has garnered much attention recently [71[8]. With 
decades of temperature measurements in hand, we know 
that this lack of correlation occurs only 0.03 — 0.1 per 
cent of the time in ACDM realizations. 

These large scales are also where cosmic variance, 
rather than statistical errors, is the limiting factor in our 
ability to compare the observed value of to its theo¬ 
retical value. This means that additional measurements 
of the temperature fluctuations will not help us make 
more definitive statements about the nature of the lack 
of correlation, and whether it is a statistical fluke within 
our cosmological model or due to unknown physics. Work 
has been done recently to quantify the viability of using 
cross correlations of temperature with E-mode polariza¬ 
tion [9] and the lensing potential ip HO] to test this “fluke 
hypothesis.” 

Correlations of CMB polarization itself, outside of just 
cross correlations with the temperature observations, are 
a natural next step in determining the nature of the lack 
of temperature correlation seen at large angles. A feature 
that is required for a real-space correlation function is 


for the field to be calculated using only local operators 
on directly observed Q and U polarization maps. The 
very nature of a correlation function that has a clearly 
defined physical interpretation depends on points on the 
sky being determined independently of each other (i.e. 
locally). 


To accomplish this, we calculate two sets of polariza¬ 
tion correlation functions: Q and U auto-correlations 
along with ^(n) and B{n) auto-correlations. These have 
a number of properties that make them unique tests of 
large-angle correlation suppression, such as contributions 
from the reionization bump that appear in polarization 
power spectra at < 10 that dominate the large-angle 
Q and U functions. The local E- and B-mode correla¬ 
tions are instead dominated by large multipoles at large 
angles, and have small contributions from reionization 
which makes them a cleaner test of physics at the last 
scattering surface. In this work we present the local 
C^^{0) and C^^(<9), along with C^^{0) and C^^(<9), 
and show distributions for the corresponding S 1/2 statis¬ 
tic for each. These results are drawn using constrained 
temperature realizations, meaning they are consistent 
with the observed power spectrum within instrumental 
errors and have a cut-sky S 1/2 at least as small as our 
cut-sky measurement. 


This paper is organized as follows: in Section [ll] we 
present the theoretical background for C{0) and a com¬ 
monly discussed statistic 5'i/2, in Section IIIII we discuss 
our calculation of the error based on next-generation 
satellite specifications as well as the lowest possible ex¬ 
pected instrument-limited value of 5'i/2, in Section ! 


IV 


present the local E- and B-mode correlation functions, in 
Section |V| we show auto-correlation functions for Q and 
U Stokes parameters, and in Section |Vl] we present our 
conclusions and discuss possibilities for future work. 
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II. BACKGROUND 


B. Stokes Q and U Correlation Functions and 
Statistics 


A. Temperature Correlation Function and 
Statistics 


The information contained in CMB temperature fluc¬ 
tuations is often represented in harmonic space by decom¬ 
posing them in terms of spherical harmonics and their 
coefficients, 


AT(n) 

To 


0(n) = ^aJ^Yem{n), 

£,m 


( 1 ) 


with the temperature power spectrum being constructed 
from the airn coefficients: 

{dim^i'm') ~ (2) 


In real space, the CMB temperature fluctuations, 
0(n), can be represented as a two-point correlation func¬ 
tion averaged over the sky at different angular separa¬ 
tions: 


Linear polarization is typically described by two quan¬ 
tities: the Q and U Stokes parameters in real space, 
and E-modes and B-modes in harmonic space. In 
real space, C^^{0) = (Qr(ni)Qr(n 2 )) and (0) = 
{Ur{ni)Ur{n 2 )) are the Q and U correlation functions, 
where Qr(n) and I/r(h) are the Stokes parameters de¬ 
fined with respect to the great arc connecting rii and 
n 2 US- Q(n) and t/(h) fields on the sphere are defined 
such that they are connected by a great arc of constant 
(j). In practice, the correlation functions are calculated as 
an average over pixels separated by an angle 0: 

C^^iO) = Qr{ni)Qr{n 2 ), 

C^^id) = Urini)Ur{n 2 ). (7) 

The decomposition of polarization into spin-2 spherical 
harmonics is done with a linear combination of the Stokes 
parameters, 

{Q{n) ± iU (fi)) = ^ ±2af„ ±2Ye^(n). ( 8 ) 

£m 


C^^{0) = 0(ni)0(n2) with hi • n 2 = cosO. (3) 

This is an estimator of the quantity = 

( 0 (hi) 0 (h 2 )), where the angle brackets represent an en¬ 
semble average. The sky average over the angular sepa¬ 
ration can be expanded in a Legendre series, 

Cf^PeicosO), (4) 

where the Ci on the right-hand side of Eq. are the 
pseudo-C^ temperature power spectrum values. 

The 5 'i/ 2 statistic was defined by the WMAP team to 
quantify the lack of angular correlation seen in tempera¬ 
ture maps [6]: 

^1/2 = J ^ d(cos6>)[C^^(6>)]2. (5) 

The expression for 5'i/2 can be written conveniently in 
terms of the temperature power spectrum and a coupling 
matrix 

si^2 = E criu'cr. (6) 

£=2 

A full expression of the matrix can be found in Ap¬ 
pendix B of [8]. The Ci fall sharply and higher order 
modes have a negligable contribution to the statistic, so 
choice of an appropriately large value of in Eq. § 
will ensure that the result is not affected by including 
additional higher-^ terms. 


The standard E- and B-mode coefficients are combina¬ 
tions of the spin-2 harmonic coefficients, 

^£m 2 {.‘^^£m —‘^^£m\ 

^£m = [2(^£m + -‘^^£m] ’ 

and the E- and B-mode power spectra are defined as 

{^£m^ £'m') ~ ^££' 

=^££'^mm'Cf^. (10) 

Using these equations, we can construct and 

from and [U: 



[Cf^Gj^icosO) -h Cf(cos6>)] , (11) 

where 

GUcosd) = - + E^p;n(cos0)) 

POS f) 

+ (i + ni)—Pr_,{cose), 

sm 0 

TYl 

((^- 1) COS ( 9 P /^( cos ( 9 ) 

— (^ -h m)P/!li(cos6>)) . (12) 

The G^(cos6>) are complicated functions of Legendre 
polynomials, so the calculation of and is not a 
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straightforward analog to Eq. (§. Instead, there will be 
three terms: 


Real-space fields of E- and B-modes are occasionally 
presented as spin-zero quantities m, 


qQQ _ riEE Ai)r<EE , A^)r<BB 
^1/2 - > . O,, 


1=2 




£ ^££' 


r,^EE A2)^BB 
■ ^££' 5 


(13) 


£m 

£m 


(19) 


where for the and are swapped. Eull details 
of calculating the matrices is outlined in Appendix H 


C. E- and B-mode Correlation Functions and 
Statistics 

The local correlation functions on the sky of the E- 
and B-modes are defined as 


C^^{e) = (B(ni)B(fi2)) 
C^^{e) = (E(ni)E(n2)). 


(14) 


The E{n) and B{n) functions can be calculated from the 
observable Q and U fields using local spin raising and 
lowering operators 8 and 3 m- 


B(n) = — [d^(Q(n) + i[/(n)) - d^(Q(n) - i[/(n))] 

-^(n) = ^ [§^(Q(n) + if/(n)) + S^(Q(n)-if/(n))] , 

( 15 ) 


where 


= — (sin6>) 

= -(sin(9)-^ 


9^ ^ \ sin 0 ) d(j) 


9 


d 


^sin^ J d(j) 

in real space, and in harmonic space. 


(sin6>) 
(sin6>) 


3 s^£m — \/ {^ — S){1 S 1) 

3 = —V"+ S){1 — S -h 1) s-lY£m 


(16) 


(17) 


In terms of spherical harmonics and coefficients, E{n) 
and B{n) are pTl 113] : 


^(“) = E o-emYemin) 

£m V E 

= E afmYemin). (18) 

£m V ' 

The prefactor under the square root is proportional to 
and is a direct consequence of using the local operators 
on the Q and U maps. 


The fields in Eq. (19) cannot be constructed from real- 


space maps only, unlike Eq. (19), and require map filter¬ 
ing in harmonic space to separate the E- and B-modes. 
Because polarization is inherently a spin-2 quantity and 
an integral over the full sky is required to extract the 
af^ and af^ coefficients from Eqs.|^and[^ the E{n) and 
B{u) are non-local. The non-local definitions of E{n) and 
B{n) require information from the full sky to separate 
the E- and B- modes from observed Q and U polariza¬ 
tion maps in any given pixel. Eor this reason, non-local 
definitions cannot be used when talking about real-space 
correlation functions, since the physical interpretation of 
a correlation at one particular point on the sky fii with 
another particular point on the sky h 2 becomes ambigu¬ 
ous. 

The expression for the two point function in terms of 
the local fields is 


= 


2£aI f{£E2)\ 
47r V(^-2)! 


Cf^P£{cosO), ( 20 ) 


and the same for the local B correlation when substuting 
in Cf^. This form of the correlation function leads to 
some interesting conclusions, namely that the traditional 
mode of thinking that 0 j is not applicable. This in¬ 
tuition was due directly to the fact that falls off as 
1/^^ and the prefactor in the sum for the TT correlation 
function in Eq. © only scales like leaving the sum 
dominated by terms less than an ^max = 30. This does 
not hold for correlation functions of the E(n) and 5(n) 
functions defined in Eq. (18), and it should be clear that 


higher £ modes will contribute to the large-angle piece 
of the correlation functions. This feature was also dis¬ 
cussed in mi, where they were focused on small-angle 
correlation functions of local E- and B-modes. 

The expressions for and are similar to 

Eq. 


qXX 

^1/2 


£max ^£max 

E 


£,B=2 


{{A±A1\ 

\{£-2)\J 


^XX T 

J-£B 


({l+n\ 

\{£'-2)\) 


riXX 
O 0! 


( 21 ) 


We have chosen to calculate the 5'i/2 statistic, rather 
than generalizing to a statistic at another angle, because 
effects that contribute to polarization inside the surface 
of last scattering (namely reionization) are at a suffi¬ 
ciently high redshift that they do not significantly change 
the relevant angle where suppression is expected to ap¬ 
pear. 
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III. ERROR LIMITS ON MEASURING A 
SUPPRESSED C(e) FOR FUTURE CMB 
POLARIZATION EXPERIMENTS 


The error in Ci for a next-generation full-sky CMB 
satellite can be determined using the relation 


AC, 


2^ + 1 



dTT j 


( 22 ) 


where a is the pixel error estimate in jaK — arcmin m- 
Values for the pixel error estimates for future surveys are 
shown in Table [T| [TBHISl . 


Experiment 

ap [jjiK arcmin] 

^FWHM [arcmin] 

Planck 

120 

5 

PIXIE 

3.78 

54 

PRISM 

3.4 

2 


TABLE I: Polarization sensitivities that reflect the actual 
Planck sensitivity in CMB channels, and the design sensi¬ 
tivity for two satellite proposals. 


To find the corresponding error band in C(^), we cre¬ 
ate 10^ realizations of the spectrum assuming chi- 

squared distribution with variance including instrumen¬ 
tal error based on the values in Table HI Constrained real¬ 
izations of are generated by drawing coefficients 
using instrument noise and assuming they are coupled to 
constrained realizations of aj^. 

The constrained temperature harmonic coefficients are 
drawn such that they produce S 1/2 values that are con¬ 
sistent with calculations from data and have a spectrum 
which matches observations (the full procedure for mak¬ 
ing constrained realizations is outlined in El). The errors 
to the mean correlation function values are determined 
based on the 68% confidence levels (C.L) for the realiza¬ 
tions. Cosmic variance dominates the error bars on the 
E- and B-mode power spectra through the reionization 
bump < 10) and instrumental error from beam size 
dominates around ^ 45 for r = 0.1. 

The instrumental error enforces a limit on the smallest 
possible value for the expectation (5'i/2), even if the cor¬ 
relation function is completely suppressed. If we assume 


that the correlation functions defined in Eqs. and 20 


are noise-free and identically zero above 60 degrees, then 
the corresponding sums over the power spectra and their 
coefficients must be zero for all P,(cos^ < 1/2). For both 
sets of correlation functions, this makes S 1/2 for Q, U, E 
or B 


Sl/2 = 



In real-space, for Q 


sc^Q{d) = 


ap 




--Qv 


(24) 


pairs 


where Vpairs is the number of pixel pairs separated by 0 
and Qrms is the root mean square value of the field. The 
integral is trivial since the only 0 dependence appears in 
the expression for Vpairs- 


N, 


pairs 




The zero true-sky value of S' 1/2 is 


^rms 


cQQ ^ 3 aj, Qi 

2A^V27ri/2 ’ 


(25) 

(26) 
rms SUbsti- 


This result is the same for the U field, with 
tuted for Qrms- 

For the E-mode statistics, it is easier to calculate 6C{0) 
in .^-space: 


6 C^^{ 0 ) = 


1 


Y^SttVj 


pairs 


\ 




(27) 


This leads to 


qEE 

^ 1/2 


8(iVpi,7r)3/2 


^^^^(2^+1)) 


w 


(^ + 2)! 

(^- 2 )! 


V^f^(2P + l), (28) 


with the same result for B when Cf^ is substituted for 
and using . 

In the near term, Planck will weigh in with its upcom¬ 
ing release of polarization data. We do not yet know the 
exact noise spectra for their EE and BB observations, 
but we can make an estimate of the expected 81/2 values 
assuming cTpoi = V2 cft and using cft = 85/iK — arcmin 
from [16]. Table |I] out lines error estimates used for Planck 
in addition to PIXIE |T7| and PRISM [18], and Table [TT] 
presents all values of the Si 12 statistic that results from 
assuming there is zero true correlation at the last scatter¬ 
ing surface for each experiment. These values show that, 
when compared to the ACDM prediction of 5'i/2, pixel 
noise is not a significant source of error to quantifying 
suppression to the correlation functions in polarization. 
Systematic errors may bias measurements of 5'i/2, but 
we will not consider these here as any unresolved system¬ 
atic would only serve to increase the value of 5'i/2- Cur¬ 
rently, no full-sky polarization maps are reliable enough 
to measure the large-angle polarization functions com¬ 
puted here. 


dcos^. 


(23) 
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Experiment 

QQ/UVluK-^] 

EEluK-^] 


Planck 

1.75 X 10"® 

0.314 

0.013 

PIXIE 

1.73 X 10"® 

3.10 X 10“"^ 1.31 X 10"® 

PRISM 

1.40 X lO"*^ 

2.51 X 10“"^ 1.06 X 10"® 


TABLE II: Expected values of S 1/2 statistic from a toy-model 
map with pixel noise using sensitivites from Table [I] and as¬ 
suming complete suppression of the true correlation function 
for Q, U, E, B. These estimates account for sensitivities for 
future CMB polarization satellites. 


IV. LOCAL B{n) AND E(n) CORRELATION 
FUNCTIONS 



In order to present a meaningful correlation function 
and related statistics, we smooth the E- and B-mode 
power spectrum with a a = 2.7° Gaussian beam (which 
corresponds to a 0.02 radian beam). There are two bene¬ 
fits to this approach: it suppresses the and for 
^ > 50 which ensures that the sum in Eq. (20) converges, 
and it suppresses all pieces of the power spectrum that 
have contributions from lensing. The former is necessary, 
since even for E- and B-mode power spectra with perfect 
de-lensing, the sum in Eq. (20) doesn’t converge through 
^max = 1500. The latter is especially important since 
we wish to make statements about correlations of pri¬ 
mordial E- and B-modes. Without smoothing we would 
need to de-lens all maps before calculating statistics. At 
the smoothing level used for analysis here, lensing does 
not contribute to the calculated S 1/2 distribution. There¬ 
fore all results used here have been produced from power 
spectra that do not include lensing effects. Figs. and 
show the resulting angular correlation function produced 
from the smoothed maps, and Figs. andshow the dis¬ 
tributions of Si /2 statistics from simulations with r = 0.1 
(smaller values of r will lead to an appropriate rescaling 
of the BB distribution, but will leave other results un¬ 
changed). For a ACDM cosmology, the best-fit value of 
is 1.86 X 10^ /iK^ and for is 218.3 /iK^. 


A feature of the correlation functions of ^(n) and 5(n) 
being dominated by large multipoles, even for large an¬ 
gular scales. These functions are also not sensitive to the 
physics of reionization, which make them a complimen¬ 
tary probe of correlation function suppression to the Q 
and U correlations presented in the following section. 


V. Q AND U CORRELATIONS 

The functions described in the section above may be 
undesirable in some cases, as they require taking deriva¬ 
tives of observations. The Q and U correlation functions 
do not require derivaties, and have the added benefit that 
they are entirely dominated by the reionization bump 
terms with ^ < 10, avoiding the need for map smooth¬ 
ing or concerns about contributions to the signal from 


FIG. 1: Angular correlation function of local B-modes r = 0.1 
with cTbeam = 2.7° Smoothing. The blue shaded region corre¬ 
sponds to 68% C.L. errors, which includes instrumental noise 
for a future gener atio n PIXIE-like experiment and cosmic 
variance using Eq. (22). 



FIG. 2: Angular correlation function of constrained local E- 
modes r = 0.1 with crbeam = 2.7° smoothing. The green 
shaded region corresponds to 68% C.L. errors, which includes 
instrumental noise for a future gener atio n PIXIE-like experi¬ 
ment and cosmic variance using Eq. (22). 


lensing. 

Fig. shows the QQ and UU correlation functions for 
r = 0.1 for AGDM. The shaded regions show the 68% 
G.L. error regions for a PIXIE-like experiment plus cos¬ 
mic variance calculated using Eq. (22). There are distinct 
characteristics of the QQ and UU functions, namely that 
the UU correlation is positive for a large range of angles 
while the QQ function is negative for a large range of 
angles. Physical suppression should drive both of these 
functions to zero. It could allow one to define additional 
measures of suppression of the correlation function be¬ 
yond the standard S' 1/2 statistic. 

Figs.and [^show the S 1/2 distributions for both the 
QQ anoUU correlation functions. The ACDM value is 
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FIG. 3: Si /2 statistic distribution for the angular correla¬ 
tion function of E-modes r = 0.1 with (Jbeam = 2.7° radian 
smoothing. The blue dashed line marks the ACDM prediction 
for the ensemble average. 



FIG. 4: Si /2 statistic distribution for the angular correla¬ 
tion function of B-modes r = 0.1 with abeam = 2.7° radian 
smoothing. The blue dashed line marks the AGDM prediction 
for the ensemble average. 


shown with the blue dashed line. The expected ACDM 
value for 5'^^ is 0.0116 and for is 0.0129 fiK^. 

In order to calculate 5'i/2, the standard efficient meth¬ 
ods defined in [9] cannot be used. Typically, Eq. (§ 
is expanded to instead be a function of the Qs and a 
coupling matrix using Eq. ( [2Q| ) rather than calculating 
the integral of the square of C{0) directly. Now, since 
Eq. (11) is in terms of Gf{cos0) rather than P^(cos6>) 
as in Eq. (20), the expressions for and become 
more complicated. Appendix [A| describes a method that 
can be used to make the calculation more efficient by 
writing G^(cos^) as functions of Wigner d-matrices. 

The large-angle Q and U correlation functions being 
dominated by the reionization era, which is entirely in¬ 
side the last scattering surface, give us a window into 
the nature of temperature suppression. The large-angle 
temperature correlation function has contributions from 


FIG. 5: Angular correlation function of Q and U polarizations 
with r = 0.1. The shaded regions correspond to the 68% G.L. 
errors. The ranges include instrumental noise for a future 
generation PIXIE-like experiment and cosmic variance using 
Eq. ([22|. 


the last scattering surface via the Sachs-Wolfe effect, and 
along the line of sight via the integrated Sachs-Wolfe ef¬ 
fect. The suppression of if caused by physics 

rather than a statistical fluke, could be due to features 
localized on the last scattering surface alone or could in¬ 
clude contributions from its interior. If features inside 
the last scattering surface are suppressed, meaning sup¬ 
pression is a three-dimensional effect, this will manifest 
as suppression in the Q and U correlation functions. 

We have chosen to calculate the standard S 1/2 statis¬ 
tic, rather than generalizing to statistics at another angle, 
S{x)^ as defined in [9], since the reionization contribu¬ 
tion is predominantly at 2 : = 10, which is near enough to 
the surface of last scattering that the angular scale that 
features subtend are nearly that of those at z = 1100. 
Contributions from late-time reionization around 2 ; = 1, 
which would skew the relevant angular scale, are subdom¬ 
inant since the amplitude of the polarization signal after 
Zreion f^lls off like a“^. This leads to an overall drop-off 
in the correlation function of meaning nearby effects 
are 100 times smaller than those at ^ = 10. 


VI. CONCLUSIONS 


To address the lack of correlation in the temperature 
power spectrum at large angles in particular, we need to 
move beyond temperature data alone. We show two vi¬ 
able methods for calculating correlation functions on the 
sky that arise from polarization and presented the distri¬ 
butions for the corresponding statistics using constrained 
realizations for the E-mode contributions and the best-fit 
ACDM framework for B-mode realizations. A suppres¬ 
sion in the primordial tensor or scalar fluctuations will 
affect the features of the two-point correlation function. 
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FIG. 6: Si /2 distribution for C^^{0) with r = 0.1. The 
blue dashed line shows the ACDM prediction for the ensemble 
average. 



FIG. 7: Si /2 distribution for C^^(O) with r — 0.1. The 
blue dashed line shows the AGDM prediction for the ensemble 
average. 


meaning , local C^^{0) and C^^{0) as well as C^^{0) 
and (0)^ and their related statistical measures. This 
would lend considerable weight to the argument that the 
lack of correlation seen in is due to primordial 

physics, and is not just an anomalous statistical fluctua¬ 
tion of ACDM. 

We presented the distribution for an S 1/2 statistic for 

a C^^{0) from ACDM cosmology with r = 0.1. If future 
limits on the value of r are found to be significantly below 
this value, the results for C^^{0) will scale appropriately, 
wheras results for all other correlation functions will re¬ 
main unchanged. For C^^{0), and C^^{0), 

we considered constrained realizations, where coef¬ 
ficients were related to aj^ coefficients that match our 
power spectrum measurements and give values of 
at least as small as we observe on the full- and cut-sky. 
We showed that for a ACDM cosmology, the expected 
values of the statistics for Stokes parameter correlation 


functions are = 0.0116 /iK^ and = 0.0129 /iK^, 
and the local E- and B-mode expected values are = 

1.85 X 10^ /iK^ and = 218.3 We chose to keep 
the previously defined S 1/2 for analysis here, rather than 
generalizing to other angles than cos 60° = 1/2, as the 
dominant secondary effect on polarization signals from 
epoch of reionzation is sufficiently close to the surface of 
last scattering to not change the relevant angle of sup¬ 
pression significantly. Late-time reionization contributes 
to the signal at a level 100 times smaller than the effect 
of reionization at ^ = 10, so while those would skew the 
relevant angular scales, they are sub dominant. 

Using a polarization error estimates for Planck, PIXIE 
and PRISM outlined in Table [Ij we calculated the re¬ 
sulting Si /2 statistics from a sky with exact suppression 
above 60°. These values are presented in Table [n| We 
note that these levels are well below the ACDM predic¬ 
tions for all of the polarization correlation functions pre¬ 
sented here, and pixel noise for future experiments will 
not be a significant source of error in identifying suppres¬ 
sion. Measurement of large-angle polarization correla¬ 
tion functions will have errors dominated by systematics 
rather than map pixel noise for the foreseeable future. 

Beyond being able to confirm that the suppression of 
temperature fluctuations is unlikely to be a statistical 
fluke, polarization correlation functions will add impor¬ 
tant new information. Because the local E and B corre¬ 
lation functions are dominated by large £ values, a sup¬ 
pression in all four correlation functions would strongly 
indicate that the suppression manifests itself physically 
in real-space at large angles. The E and B correlations 
give insight about suppression that is independent of any 
effects of reionization which dominate the Q and U cor¬ 
relations. Also, foreground emission will contribute dif¬ 
ferently to the various correlation functions. 

Eurther, since the local B correlation is determined 
entirely by tensor fluctuations, a strong suppression in 
that correlation function and not in others would show 
that the primordial suppression is predominantly in the 
tensor perturbations, while suppressions in local E, Q 
and U but not in local B would suggest that the scalar 
perturbations are suppressed. 

The distribution for S 1/2 statistics for each constrained 
correlation function was compared to the distribution 
from ACDM alone. We found no significant difference 
between the two distributions and have presented only 
the constrained in this work. This means that polariza¬ 
tion correlation functions provide a largely independent 
probe of correlations compared to the anomalous temper¬ 
ature correlation function. Euture high-sensitivity mea¬ 
surements of polarization over large fractions of the sky 
from envisioned experiments like PIXIE m will differ¬ 
entiate primordial physics from a statistical fluke as the 
origin of this anomaly. 

If the suppressed temperature correlation is due to a 
statistical fluke, then measurements of the polarization 
correlation function at large angular scales is likely to 









give a much less suppressed sigual. If, ou the other haud, 
the suppressed temperature correlatiou is due to some 
physical mechauism, how well cau polarizatiou test this 
sceuario? The auswer depeuds ou the precise predictiou 
of the suppressiou model. A Bayesiau model comparisou 
betweeu a giveu model aud the staudard cosmology will 
give a quautitative auswer to this questiou. We are cur- 
reutly iuvestigatiug this possibility for a model with sup¬ 
pressed correlatious iu the primordial gravitatioual po- 
teutial perturbatious. Iu geueral, we expect suppressed 
primordial correlatious will be evideut iu polarizatiou 
at least as much as iu temperature, due to the lack of 
au iutegrated Sachs-Wolfe coutributiou to the polariza¬ 
tiou perturbatious. A stroug discrimiuatiou betweeu a 
suppressed-correlatiou cosmology aud the staudard cos¬ 
mology is likely. 
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Appendix A: Correlation Functions and S(x) 
Calculations for QQ and UU in Terms of Wigner d 
Matrices 


11 


cau be ex- 


The Gf {cos 0) fuuctious described iu 
pressed iu terms of reduced Wiguer matrices. This form 
may be useful for fiudiug aualytic expressious of 5'^^ aud 


^ 1/2 which are easier to calculate uumerically thau per- 
formiug the full iutegrals over [C{0)]‘^. We cau write the 
correlatiou fuuctious for Q aud U as 


CQQ^O) = ^ {D+{cosO)Cf^ + D^icosd)Cf^) 

£ ^ 

(Al) 


aud 

^ {D^{cose)cr + Di{cose)CD , 

£ 

(A2) 

where we have assumed parity iuvariauce aud used 

DficosO) = [dl^ie) ± 4,-2 (^)] (A3) 

from m, where d\ -{0) are the reduced Wiguer rotatiou 
matrices. 

This form of the correlatiou fuuctious would lead to a 
method of calculatiug Si 12 most similar to that defiued 


iu [9] by usiug properties of d-matrix iutegrals aud recur- 
siou relatious. 

The geueral form of the S statistic is defiued as 

S{x) = J [C^^{x)f Ax, (A4) 

where x = cosO aud Si 12 = 5'(l/2). We defiue 

= J 4.±2(a:)4U2(a;) da;, (A5) 


which we use to calculate S{x) usiug properties of the 
reduced Wiguer matrices. 

Importaut properties of the matrices are that 

^u~{x) = l^^{x) aud 


r,,7{x) = (-1)^+^' 


2 

2 ^ + 1 


kv - 


, (A6) 


SO all required quautities cau be coustructed from calcu¬ 
latiug /^r (x) aud (x) ouly. To compute these, we 
use the relatiou betweu reduced Wiguer matrices aud the 
Clebsch-Gordau coefiicieuts [20], : 


,1712 ,7712 (^) 

£+£' 

j=\£-£'\ 


Combiuiug Eq. (A5) with Eq. (A7) aud exploitiug prop¬ 
erties of the Clebsch-Gordau coefficieuts, we hud: 



E 


jf=max(|^—1,4) 



d{^{x) dx (A8) 


aud 


(a;) 


£+£' 

E 


j=ma.:x.{\£—£' \ ,4) 



4o(a:) da;. 


(A9) 


The iutegrals have aualytic solutious, which make com- 
putatiou of the l££f matrices more efhcieut. The iutegral 
over d 4 Q(x) is the simpler of the two cases: 



(AlO) 


This iutegral cau be performed by uotiug that 


dLM = 

(All) 

Usiug the Rodrigues formula, 

Pf(a;) = (-l)-(l - x^r/^—P^{x), (A12) 
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and integrating by parts, we can show that 


where 



+ AxPf{x) - 4(1 - 

ax 

P 4.{-iy j{j P 1) - 24xPj{x) - 24{-iy 
24 \ 

+ 27Ti • (A13) 


kj{x) = 


£j{x) = 


l{j + 2){j-l) 
(j-2)(i + 3)^^ 
2 




H- c?42(a:), (A17) 

V(i-2)(i + 3) 

4o(a;) - 


VjU + 1) 


\/i(j +1) 


mj{x), 


and 


(A18) 


The integral over P^^{x) is more conveniently done as 
integrals over angles: 

if\x)= j d^ 4 (^)sin^d/ 3 . (A14) 

J cos“i (a:) 

Using the relation 


sin p 


^(j2-42)(j- + 3)(j+4) 

j(2i + i) m 


4v/(i-3)(j + 4) 

j(i +1) 




^[(j + 1)2-42](j--3)(j + 2) 

(j + l)(2j + l) "3 A). 

(A15) 


along with properties of integrals over 6 / 43 (;d), d\ 2 {P)i 
and d\Q{y)/ sin(/3), the integral becomes 


mj{x) = 


+ 


/ (j-4)! 

/(j + 4)! 

2 x 


^/l — x‘^ 


Pf{x) 


1 — X 


Pix) + 


2 3 


(-ir (j + 2)! 

4 (j - 2)! 


-2—P,(:r)-(-iri(j + l) 


(A19) 


Recursion relations were used to calculate the Clebsch- 
Gordan coefficients from Eqs. (8.5:3), (8.5:8), and 
(8.6:27) in [20]. 

We compute the matrices in Eq. ( [l3| directly from 
these forms, as 




j{2j + 1) 

4y(j-3)(i + 4) 


kj-i{x) 


kj{x) 


jU +1) 

^[(j + 1)2-42](j-3)(j+2) 

(j + l)(2j + 1) 


kj+i{x), 

(A16) 


I^l}{x) = I++ix) + 21+7 (x) + Iu7{x) 

lf^{x) = I++{x)+Iu7{x) 

ke' = Ip' (x) - 21 jf {x) + {x). (A20) 
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